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We present a model-independent, non-perturbative relation between finite-volume matrix elements 
and infinite-volume 0 —>■ 2 and 1 —>■ 2 transition amplitudes. Our result accommodates theories 
in which the final two-particle state is coupled to any number of other two-body channels, with all 
angular momentum states included. The derivation uses generic, fully relativistic field theory, and 
is exact up to exponentially suppressed corrections in the lightest particle mass times the box size. 
This work distinguishes itself from previous studies by accommodating particles with any intrinsic 
spin. To illustrate the utility of our general result, we discuss how it can be implemented for studies 
of A -t- y — >■ [Ntt , Nrj, Nff ,Y1K, KK) transitions, where A is a generic external current. The 
reduction of rotational symmetry, due to the cubic finite volume, manifests in this example through 
the mixing of S- and P-waves when the system has nonzero total momentum. 
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I. INTRODUCTION 

Whether catalyzing the formation of exotic mesons (as is the case for the GluEx experiment mm) or probing 
the Standard Model of particle physics (e.g., heavy meson decays [SHHl), hadronic transitions lie at the core of 
present-day nuclear and particle physics research. Consequently, there is great incentive to provide theoretical pre- 
and postdictions of these reactions. However, these quantities are very challenging due to the nonperturbative nature 
of the strong nuclear force and its underlying quantum theory, quantum chromodynamics (QCD). The preferred 
approach to overcome this obstacle is to implement numerical studies of QCD, in particular lattice QCD (LQCD). 

Numerical LQCD calculations are necessarily performed in a finite Euclidean spacetime. For processes involving 
two or more hadrons in the initial or final states, the relationship between quantities determined from LQCD and 
the infinite Minkowski spacetime observables is therefore highly nontrivial.^ Nonetheless much progress has been 
made by using finite volume as a tool, rather than an unwanted artifact, in order to extract physical observables. In 
particular, the problem of extracting transition amplitudes from finite-volume matrix elements was first addressed 
by Lellouch and Liischer in the context of A —>■ tttt decays [2D]. This idea has since been revisited and generalized 
to describe more complicated one-to-two processes jUMj, two-to-two processes [251430] . as well as non-relativistic, 
weakly interaction multi-bosonic systems |31| . Related ideas have also been presented to describe the finite-volume 
effects for the long-distance two-pion contribution to neutral kaon mixing [32] . 

In this work, we extend the ideas presented in our previous article. Ref. HZ], and give a model-independent relation 
between finite-volume matrix elements and transition amplitudes for particles with arbitrary intrinsic spin. We 
primarily consider matrix elements of local currents coupling one- and two-particle finite-volume states and relate 
these to infinite-volume 1 —>■ 2 transition amplitudes in the presence of an external current. We additionally consider 
finite-volume matrix elements connecting the vacuum to two-particle states and relate these to infinite-volume 0 —> 2 
transitions. Here we also present an alternative derivation from our earlier work. Instead of considering ratios of 
two- and three-point correlators to access the desired matrix elements, here we access them directly from a single 
finite-volume two-point correlator. This is possible because we have sufficient freedom for the operators appearing 
within the correlator and can choose these to define both 1—^2 and 0 —>■ 2 transitions. 

As was the case in Ref. ^ 7 \, our result here allows for two-particle states that strongly couple to any number of 
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open two-particle channels and also incorporates all angular momentum states. The only approximation required is 
neglect of exponentially suppressed corrections of the form where m is the physical mass of the lightest particle 

and L is the extent of the finite volume. Therefore our results apply only to theories where mL ^ 1 and do not, for 
example, describe QED-I-QCD in finite-volume.^ 

Beyond facilitating the determination of a wide range of physical observables, this work will also impact the 
assessment of systematic error for quantities that have already been extracted from LQCD calculations. One prominent 
example is the study of parity violation in the nuclear sector performed by Ref. [39j . This work constrains “ 

1.099(51)(6) X 10“^ via LQCD using light quark masses corresponding to ~ 389 MeV. This is the coupling that 
parametrizes the parity violating N —>■ ttN process at threshold, where the final two-body state is in an S-wave. At 
the time of this benchmark calculation, it was not understood how to correct finite-volume effects due to the fact that 
the final state is composed of two particles. In Section |V] we explicitly discuss the implication of the main result to 
N + J ^ {Ntt, Nr], Nr]', EAT, AK), where J is any external current. 

We comment that, although this work is primarily intended for application in numerical LQCD, the result presented 
is universal and holds for any quantum theory in a finite volume. For example, Ref. |40j studied radiative neutron 
capture (n -|-p —> d + "/) in lattice effective field theory (EFT), where nucleons are treated as fundamental degrees 
of freedom. ^ By considering volumes where 3> 1, where Kd is the binding energy of the deuteron and L is the 
spatial extent of the volume, one may use the formalism presented to study this and many other reactions (e.g., 
p + p ^ d-ke"*" + Ue or d + u p + n + v). 

The remainder of this work is organized as follows. In the following section we completely specify our set-up and 
notation. We then define a finite-volume correlator and state a key identity, expressing it in terms of infinite-volume 
quantities and finite-volume kinematic functions. In Section [TTT| we review the derivation of this key identity. Then, in 
Section EYi we show how the identity for the finite-volume correlator can be used to derive all results of this work. In 
Appendices [A| and |B| we discuss two important simplifying limits of this result, namely the free limit and the narrow 
width approximation. In Appendix we explain how to generalize this result to volumes that are rectangular prisms 
with twisted boundary conditions. 

We close this section with a summary of our main results, all derived in Section [IV| using the expression for the finite- 
volume correlator, Eq. (35). We first show how the expression can be used to relate the finite-volume energy spectrum 


to two-to-two scattering amplitudes. The relation takes the form of a quantization condition nanaiisiiaiiiHg 


A(P, L) = det[F-^(P, L) +M{P)] = 0, 


( 1 ) 


where P is the four-momentum of the two scattering particles. The interpretation is that, for fixed total spatial 
momentum P and box size L, the energies E 2 ^p^l, • • ■ for which A vanishes are precisely the energies in the 

spectrum of the finite-volume theory. Since A is a known kinematic function [see Eq. (33)], if the energy spectrum is 
also known, then this result can be used to constrain the scattering amplitude, A4. The notation needed to completely 
understand this result is discussed in Section [Tl] In arriving at this result, we assume the two-particle states lie below 
the multi-particle thresholds. For recent progress towards studying three-particle systems in a finite volume we point 
the reader to Refs. [iKM7] . 

We next derive the relation between finite-volume matrix elements and 1 —> 2 transition amplitudes [20H231E5H27] 


I«, P', L|J)4(0, P - P') I Ao,p, P, L, 1) I = nK, P') > 

where Afi and Afn are the normalization factors for the one-body and two-body finite-volume states. 


( 2 ) 


(Eo,P,P,L,l|Ao,P,P,L,l) =A/'i, 

(£;;,p',lk,p',l)=aa„. 


(3) 

(4) 


We set these equal to one for the remainder of this work. Here we have also introduced 

Tl{En,P) = lim ' 


— {iPi + En)-^ 


F-i(P,L)+M(P) 


(5) 


^ When performing numerical lattice QCD+QED studies, one may remove the zero mode, which allows for a systematic assessment of 
finite-volume effects |33H37| including for systems involving two-particles |38| 

^ In their work, Rupak and Lee used an alternative technique to the Lellouch-Liischer formalism in determining the transition amplitude. 
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[nT{Eo,P: P; E'^, P')] (27r)353(P - P' - Q) = (i?;, P', { J'}, out|J^(0, Q)|i?o,p, P, {^}i) ■ (6) 

Here and below we use En in place of to reduce clutter of notation. In Eq. ^ we have introduced the labels 

{J}i and {J'}, which stand for lists of all quantum numbers specifying the one-particle in-state and two-particle 
out-state respectively. For the out-state, {J'} includes a channel label denoting the flavor of the individual outgoing 
particles. In Eq. ([^, is understood as a row, TZ as a matrix and as a column in this index space. J'a is a 
generic external local current. The matrix TZ is the generalization of the Lellouch-Liischer proportionality factor that 
enters the relation for extracting K —>■ tttt decays. Note that evaluating the limit in practice requires determining 
the first-derivative of the scattering amplitude evaluated at the energy pole. Again, see Section [IT] for a complete 
explanation of notation. 

Finally, Section |IV| contains the derivation of two additional results. First we present a method for constraining 
the relative sign of two different transition amplitudes, involving the same in- and out-states but different current 
insertions 


(K, P', L\Ja, (0, P - PQIAq.p, P, £, 1} 
(f;, P', L\JaA0, P - P')|i^o,p, P, L, 1 ) 


At 7^(F;;,p') nz 

At 7^(^;;,p') nz 


(7) 


where At is a completely arbitrary vector to be chosen at the user’s convenience. This result follows from the 


observation that TZ can be written as an outer product of row and column vectors [see Eq. (100)] below. This outer- 
product form suggests an interpretation of Eqs. ([^ and Q as relations between finite- and infinite-volume two-particle 
states, an idea already introduced in Ref. [20] and also discussed in Refs. [Ml 1^ . 

Second, we derive the results for 0 —>■ 2 transitions, which are reached by a straightforward modification of those 
just given. We find 


where 


|(F„,P,L|J^(0,-P)|0)| = 



L37^(£;„,P) 


V 


out 
A ’ 


[V1"‘(F„, P)] (27r)353(p ^ Q) ^ out|JA(0, Q)|0). 


( 8 ) 

(9) 


We emphasize that the formalism presented here only holds for local currents and cannot be implemented for processes 
that contain long range contributions, e.g., 7 * 7 * —?> tttt. 

This equation will allow for future studies of decay constants of hadronic resonances. One example that has received 
attention in the lattice QCD community is the decay constant of the 7r(1300) resonance jlSJUnj. To be able to study 
this resonance at or near the physical point would require solving the many-body, multi-coupled-channel problem in 
a finite volume. This is beyond the scope of current calculations due to the issues discussed above. For unphysically 
heavy light quark masses, however, the possible strong decay channels reduce dramatically. For instance, at the SU(3) 
flavor point, studied in Ref. [15], the only open channel is the p 7 r-scattering channel, where the p is itself stable under 
the strong interactions. In this limit, one can utilize Eq. Q to determine the p 7 r —>■ 0 transition amplitude. From 
this, one may analytically continue the 7r(1300) pole, and thereby determine its decay constant. 


II. IDENTITY FOR THE FINITE-VOLUME CORRELATOR 


In this section we review the decomposition of the two-to-two finite-volume correlator, denoted C'l(P), into products 
of finite- and infinite-volume quantities. This decomposition was presented for scalar particles in Ref. |22] and for 
the Vtt system in the A-resonance channel in Refs. EaEDI- The result reviewed here accommodates particles with 
arbitrary spin. We begin by describing the set-up of the calculation as well as the notation that we use. 


In this work we describe particles and interactions using a completely general relativistic field theory. [See Eq. (36) 
below and the discussion that follows.] To restrict attention to two-particle states we require that the center-of- 
mass (CM) frame energy is below the lowest threshold with more than two particles. Thus, all on-shell states contain 
exactly two particles. Here we accommodate N,. different channels each containing a particle pair with given individual 
physical masses and spins. We use the index a = 1, • • • , A), to label the channels considered; TOqi, 771^2 denote the 
two physical (pole) masses and s“, S 2 denote the two spins. 

This theory is then considered in a finite-volume. In the present work we restrict attention to the simplest possible 
implementation, a finite cubic spatial volume with linear extent, L, and with periodic boundary conditions. The 
generalization to asymmetric boxes and to twisted boundary conditions is well understood and we give the equations 
in Appendix We further assume that the time extent is taken to be infinite, and that L is large enough such that 
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exponentially suppressed corrections of the form 6“"*^ can be neglected, with m the physical mass of the lightest 
particle in the spectrum. In our analysis we allow for nonzero total momentum in the finite-volume frame. We use P to 
denote the total momentum in this frame and E to denote the total energy. Due to the periodic boundary conditions, 
P must equal a vector of integers multiplied by (27r/L) [P £ (27r/L)Z^]. Finally we introduce E* = \/E'^ — for 
the CM frame energy. As already mentioned, E* is constrained to a range of energies for which only two-particle 
states can go on-shell.^ 

Within this finite-volume set-up we define the correlator 


Cl{P) = Cl{P 4 ,P) = J^d 




(o|rM(x)si'(o)|o) 


J L 


( 10 ) 


Here T indicates time ordering and A{x), B^{x) are operators which couple to the two-particle states of interest. The 
general allowed forms for A{x) and B^x) are of great importance to the arguments made in Section IV We discuss 


the general conditions for which our analysis of Cl holds in the text following Eq. (49) below and again in Section 


|IV| where specific choices for the operators are considered. We use the Euclidean signature metric throughout with 
Xi = ix^ and P 4 = iP^ = iE. The subscript L on the correlator indicates that it is to be evaluated using the Feymann 
rules of finite-volume field theory. As was shown in Ref. [52], one can express Cl{P) in terms of the infinite-volume 
correlator Cao{P), infinite-volume matrix elements of A(0), B^{0), infinite-volume two-to-two scattering amplitudes, 
and known finite-volume kinematic functions. We carefully define the matrix elements, scattering amplitudes and 
kinematic functions and then state the result for Cl{P) in Eq. ( |3^ below. 

To define the infinite-volume matrix elements of A(0) and ^8^(0), as well as infinite-volume scattering amplitudes, 
we first need to introduce notation for two-particle states with arbitrary spin. We first consider two-scalar states. 
These may be specified by total energy E, total momentum P, momentum of a single particle k and channel a: 
|i?,P,a, k). The states are assumed to satisfy standard relativistic normalization 

{E', P', a', k'lE, P, a, k) = 5aa'2uJia2uj2a [5^(k - k')(53(P - k - P' + k') + <5(a)<53(k - P' + k')(53(P - k - k')] , (11) 
where i5(a) = 1 if the particles are identical and 0 otherwise. We have also introduced the shorthand 


UJal = \J+ k^ , UJa2 = \/™a 2 + (? “ k)^ ■ (12) 

It is important to note that, with E and P fixed, k cannot be freely chosen. Indeed the only remaining degree 
of freedom for two on-shell scalars is the direction of motion for one of the two in the CM frame. To see this, note 
that if k is chosen to satisfy the total energy condition E = uJai + ^a 2 : then a boost of the four-vectors (w^i, k) and 
(a;a 2 ,P — k) to the CM frame, i.e. with boost velocity /3 = —PjE, gives the four-vectors (\/p1oi 
(\/”^ a 2 + ~9a^a) respectively. Here g*, the magnitude of particle momentum in the CM frame, is given by 

E* = \lml^ + qf + a/ ml^ + qf . (13) 


We deduce that k^ is the only remaining degree of freedom for fixed A, P so that we can rewrite the two-scalar state 
as |E, P,a,kjj). It is further convenient to decompose the two-particle states in spherical harmonics. In this work we 
will need to consider both in- and out-states and we thus define 

|A,P,a,k*,in) = (kl)|A, P, a, / ,TO/,in), (14) 

(F;,P,a',k*,,ont| = V47rl)T™; (k*/)(E, P, a', Z',m), out|. (15) 

To include spin we simply add the labels s“, to the two-particle states. Here mg- labels the azimuthal 

component of the spin of the jth particle, with respect to some arbitrary fixed axis in that individual particle’s CM 
frame. A particularly convenient choice is to quantize the spin of the particle along its finite-volume-frame momentum. 
This defines the helicity, which we label as Ag^. Since the full set of indices is long we also introduce the shorthand 


|E,P,{;},in) = \E,P,a,l,mi,sl,Xs,,s^,Xs 2 ^E) , 
(£;,P,{Z'},out| = (A,P,a',;',m;,s'i“, A(,^,S 2 “, A(,^,out|. 


(16) 

(17) 


4 


In the case of a Z 2 symmetry that decouples even- and odd-particle-number states, the required restriction is 0 < E* < 4m. 
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An alternative basis for two particles with spin is reached by replacing the indices Z, mp s“, , S 2 , with 

J, M, l,S,Si,S 2 - Here J is total angular-momentum in the CM frame, M is the azimuthal component of total-angular 
momentum in the CM frame and S is the total spin. The two bases are related by 

|A,P,{J},in)= ^ |A,P,{Z},in)({Z}|{J}), (18) 

Ajj ,As 2 .mi ,ms 

where we have introduced analogous shorthand for the other basis 

|£;,P,{J},in) = |T;,P,a, J,M,Z,S',s)‘,s^,in), 

(A, P, { J'}, out| = {E, P, o', J', M', I', S', out|, 

and where 

({^}|{-^}) = {siXsi,S2Xs2\Sms){lmi,Sms\Jmj), ( 21 ) 


(19) 

( 20 ) 


with the right-hand side equal to a product of helicity Clebsch-Gordan coefficients. Note that the Clebsch-Gordan co¬ 
efficients are more commonly written in terms of states whose spin is quantized along the £-axis, (s“ , S 2 Ws 2 15' rns). 

One can write ({Oil'll) terms of these and Wigner-I? matrices. Let Rj be an active rotation from (0,0, m |) to k_,-, 

(s') 

where kj is the momentum of the jt/i particle. Then define 4. (i? 2 ) asthemg^Aj component of the corresponding 


Wigner-I? matrix in the Sj representation. With this, we find that Eq. (21) can be written as 


E 


(s? ms ,, 4 TOg^ \Sms){l mi,Sms\J (i?i)T>^(^4 {R 2 ) ■ 


( 22 ) 


Having completed our introduction of two-particle states, we now define the infinite-volume matrix elements of 


A(0) and H^O) which enter the identity for Cl{P) 

A{q(P) = A{q(P4,P) = (0|A(0)|-zP4,P,W,in), (23) 

p|j,j(P) = P|,,j(P4,P) = (-zP4,P,{/'},outlet(0)|0), (24) 

or in the total-angular-momentum basis 

Ap}(P) = A{j}(P 4,P) = (0|A(0)|-iP4,P,{J},in) , (25) 

P},,j(P) = P},,j(P4,P) = (-*P4,P,{J'},out|pt(0)|0). (26) 


A subtlety arises at this stage. The states appearing in the equations above can only be interpreted as two-particle 
states for —iP 4 = E real and positive. However the field theoretic definition of the states can be extended into the 
complex plane so that A(P) and B'^{P) can be defined for all values of P 4 . We assume that this extension has been 
performed but it will play little role and will not enter our final result. Finally, below we will often suppress the 
indices and write A(P) to represent a row-vector and B''{P) to represent a column. In all expressions where this is 
done the equation is correct in both and {J} bases. 

We turn next to the two-to-two scattering amplitudes, collectively denoted M., which enter the identity for Cl{P)- 
These are given by 

(27r)44(P-P') = (-^P(,P^K},out|-^P4,P,{0,in)eon„, (27) 

(27r)44(P-P') = (-iP.^,P^{T'},out|-^P4,P,{J},in)eo„n, (28) 

where the subscript “conn” indicates that only fully connected diagrams should be included. These scattering ampli¬ 
tudes are related to the S-matrix via 


iM(P) = P-4S'(P) - l]p-4 

Here I is the identity element, corresponding to the disconnected diagrams that we omit, and 

p = ^^^;^ diag . . . , ViNq*N) , 


(29) 


(30) 


with equal to 1/2 if the particles in the oth channel are indistinguishable and equal to 1 otherwise. For —iPj = E 
real and positive, this is the standard definition of the two-to-two scattering amplitude and S-matrix. Analytic 
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continuation to the entire complex P 4 plane is also well understood. As with A{P) and B^{P), we will often suppress 
the indices below and write M{P) which is understood as a matrix. 

Finally, the kinematic quantities which enter the identity for Cl{P) are 


F{j},{r}{P,L)^5ss' Y. {lm,Sms\JM){l' m[,S' ms\J'M')Fl^^^^,i,^,{P,L), 


mi ,m' ,7725 


where 





,(kj 


2uJal 2Wo2 ( a - Wal - a;a2 + *e) V 9a 
Here the superscript “sc” stands for scalar. We have also introduced the notation 


i+i' 


' ^ 'sf'" 

- 


. 


1 

Z3 


E 

kG(27i-/L)Z3 


dk 


(31) 

(32) 


(33) 


(34) 


In Eq. (33), k* = A:*k^ is the spatial component of (a;*^,k*), which is the four vector reached by boosting (a;ai,k) 


with a boost velocity (3 = —VjE. Note that the magnitude fc* is not equal generally equa l to g* introduced above, 
because, for a general choice of k, the condition E = Wai + Wa 2 is not satisfied. See Eq. (1241 below for the alternative 
definition of that was first given in Ref. [22] . 

Our claim, embodied in Eq. is that spin can be incorporated by trivially modifying the results for scalar 

particles derived in Ref. jH]. We demonstrate in the next section that this straightforward modification, simply 
including Kronecker deltas in helicity, is the correct prescription. Using the total angular momentum basis, Eq. (32), 


results in a more signihcant difference in the form of F between scalar and non-scalar particles. The relations between 
the two basis follow from applying the transformations specihed in Eqs. (18) - (22). Note that the distinction between 


helicity quantization and spin quantization along some arbitrary axis does not enter into Eq. (32). In other words, 
the Wigner-D matrices appearing in Eq. (22) cancel out in the definition of 
We are now ready to state the key identity for Cl{P) 


Cl{P) = C^{P)-A{P) — 


1 


F-i(F,L)+M(P) 


B^P). 


(35) 


Since A{P) is a row-vector, A4{P) and F{P,L) are matrices, and B\P) is a column vector, the second term of 
Eq. (35) has no uncontracted indices. It is for this reason that the second term is basis independent, and can be 


expressed using both {Z} and {J} indices. In the following section we review the derivation of this result, and in 
Section [Tv| we use it to derive our relations for extracting 0 —)■ 2 and 1 —>■ 2 transition amplitudes. For the latter case, 
we reach the relation by choosing A and to contain both the current and the single-particle creation operators. It 
is for this reason that we must also carefully consider the types of operators A and for which Eq. (35) holds, to 
justify this non-standard construction. 


III. DERIVATION OF EQ. ( |35| ) 


In this section we review the derivation of Eq. (35) in three parts: 


1. Argue that Cl{P) is equal to an infinite diagrammatic series which can be organized into a skeleton expansion 
built from two-particle loops (with summed momenta) and infinite-volume two-to-two Bethe-Salpeter kernels 
[see Figure [^. 

2. Apply an identity which rewrites the summed two-particle loops as analogous integrated loops plus a finite- 
volume residue correction, given by F{P,L) [see Figure [^. 


3. Reorganize the terms in the series by the number of factors of F they contain, and then sum all infinite-volume 
diagrams which appear outside and between F factors. We find that the inhnite-volume diagrams outsid e F 
factors sum to A{P), B^P) and those between to M.{P). Summing the series to all orders in F gives Eq. (35) 
[see Figure [^. 


We now discuss each of the steps in some detail. 
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FIG. 1: (a) Depicted is the diagrammatic representation of the two-point correlation function in a finite volume for energies 
where only two-particle states can go on-shell. B"' and A denote the creation and annihilation operators respectively. The gray 
circles depict the sum of all diagrams that are two-particle irreducible in the (i?, P)-channel, defined in (b). The propagators 
are fully dressed and explicitly defined in terms of the one-particle irreducible (IPI) diagrams in (c). Although in the figure 
we only explicitly depict the contributions to the correlation function of a single channel, the formalism presented holds for an 
arbitrary number of open, two-body channels. See Ref. m for a diagrammatic representation with an arbitrary number of 
open channels. 


A. Skeleton expansion 


To reach the skeleton expansion we start with the general diagrammatic expansion which defines Cl{P)- This is 
determined by a set of Feynman rules which are in turn specified by the Lagrangian density of the theory 


Ns 


Nf _ 

£ = (1/2) ^ -I- ^ -h • • • -I- P(/>, 'F, • • •) ■ 


(36) 


i=i 


Note that all terms are written with Euclidean signature, for example d^d^. Here we have explicitly shown 

Ns different scalar particles, with corresponding fields ^i>i, • • • , , as well as Nf different spin-1/2 fermions, denoted 

dll, • • • ,'FjVj ■ To define the fermion kinetic terms we have introduced the shorthand $ = where Je,/! are 

Euclidean gamma matrices defined to satisfy {^E,tj.,"fE,v} = The first ellipsis in Eq. (361 stands for the kinetic 
terms of other particle types, with spin exceeding one half. All possible interactions are include in the potential 
E((/, dl, • • •), where the ellipsis stands for the dependence on all other particle types. These could include fields for 
spin-1 mesons (e.g., the J/dl and B*), stable spin-3/2 baryons (e.g. the H-baryon), etc. The interactions are assumed 
to be local and to respect Lorentz invariance but are otherwise arbitrary. We further assume that all counter terms, 
including mass and wave-function renormalization, are included in E((/, dl, • • •). We constrain these using on-shell 
renormalization, meaning that all m parameters in the Lagrangian density are physical pole masses and that the 
wave-function renormalization of all fields is equal to one. We stress here that this Lagrangian is understood to be 
in terms of the low-energy degrees of freedom of the theory, so that there exists a correspondence between the fields 
shown and the low-lying particle spectrum. 

The resulting Eeynman rules include standard propagators of the form 


^.free(£) “ p2 + ^2^^ ’ “ ^2 ^ 2 ^^ 


(37) 


where, as above, the ellipsis stands for the analogous propagators of higher spin fields. We will return to the specific 
forms of the higher spin propagators below. Interactions and counterterms are encoded by a set of vertices, whose 
specific form is determined by E( 0 , d>, • • •). 

To calculate C'l(P), one also requires the Feynman rules induced by the interpolators A and . These are 
determined by first expressing the operators in terms of the elementary fields appearing in Eq. (36). For scalar 
particles, for example, the decompositions generically take the form 




(38) 


and similar for . Here the ellipsis stands for terms constructed from three or more single particle fields. In each 
term the operators are integrated over weight functions, denoted Note that the weight functions cannot depend 










on X as this would violate the translation property 


g-*p.y+i^y4^(2.)g^P•y-i^y4 ^ _4(3. ^ 


(39) 


which must hold if C'l(P) is to be projected to a single total momentum. The indices ij ■ ■ ■ are applied to the weight 
functions since these can be different for different particle flavors. These are also applied to the symmetry factors 
which equal prodcuts of 1/n! factors for each subset of n identical scalars. 

Turning next to spin-one-half particles, we generalize the operator decomposition by replacing the scalar fields with 
Dirac helds and then additionally including Dirac indices on the weight functions 


-YlJ+ y) + Jd^y J + + + 


(40) 


This form readily generalizes to also accommodate higher spin particles. For a particle with spin s there must exist 
some “Dirac-like” indices which encode the spin degrees of freedom. These indices generally do not correspond to 
dehnite spin states. In the following section we will describe a change of basis that converts such field indices, which 
label the interpolating functions and propagators, to spin indices that label physical states. However, as a first step 
we express the operators in terms of indexed weight functions as in Eq. ( |40[ ). We finally note that, for our fully general 
theory, the operators A and will be a sum over all of the different spin types included in the theory. 

Having discussed the general forms of the interpolators A and in terms of single-particle fields, we are now in 
position to give the Fey nma n rul es a ssociated with these operators. These are given by Fourier transforming the 
weight functions of Eqs. (381 and (40). For example the single field terms of A and B^ in the scalar sector, first term 


in Eq. (38), induce single-legged vertices with vertex factors 


if = Jd%APyAf\y), 


d^ye- 




(41) 

(42) 


where here and below we suppress the dependence on total momentum P. Similarly, the two-field terms induce two 
legged vertices with factors 


(fc) = Jd"yJ{y, z ), 

^ J^iyJ ^(y,z). 


(43) 

(44) 


This pattern continues to all orders with the nth vertex factor equal to n Fourier transforms with momenta coordinates 
constrained to sum to P. Finally the generalization to spin is achieved by applying field indices to both sides of 
Eqs. (|4Tl)-(|44l). 


With our general Feynman rules in hand, we next envision enumerating the full set of diagrams contributing to 
Cl{P)- We then organize this series into a simplified skeleton expansion. To do so one must identify all parts of 
diagrams with two-particles that can simultaneously go on-shell. Such particles must carry the total energy and 
momentum E, P and must also have the quantum numbers of the operator B^. Suppose that a given diagram has N 
distinct pairs of propagators which carry E,P. Then this diagram is contained within the A^th term of the following 
skeleton expansion [see also Figure 


= ZuE / (fc)‘5a(fc)<f fc) 

k 


(45) 


where 


Ml 


'Ak', k) = -Ka'Ak', E J '^^a',b{k',p)Sb{p)ML,bAP’ k) . 


(46) 


^ A technical assumption first used here is that correlators and amplitudes, in both finite- and infinite-volume, are given by summing 
perturbation theory to all orders. 
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Here 


Sa{k) = ^aK\ik)KliP - k) , 


(47) 


is the product of the two-propagators in the ath channel. For particles with nonzero spin, this should be understood 
as a direct product of matrices that reside in two distinct spaces. The propagators here are fully dressed and are 
defined in terms of the two-point function as 


A:\ik) ^ , 


(48) 


where is an interpolating held for a spin s particle and is the appropriate conjugate held. For example, for 

complex scalar helds and for Dirac fermions 

The functions Ka'^aik', k) are two-to-two Bethe-Salpeter kernels with incoming particles in channel a with momenta 
k, P — k, and outgoing particles in a' with momenta k',P — k'. The Bethe-Salpeter kernel equals the sum of all 
diagrams that are two-particle irreducible with respect to propagator pairs carrying (i?, P). Equivalently, the kernels 
contain all diagrams with no on-shell intermediate states. Following Ref. [22] . on can prove that for loop momenta 
in diagrams with no on-shell states, the difference between sums and integrals is exponentially suppressed as 
Since we neglect these corrections, we work throughout with the inhnite-volume (integrated-loop-momenta) forms of 
the Bethe-Salpeter kernel iFa',o- We also use inhnite-volume propagators in Sa, which is justihed since self-energy 
diagrams do not contain on-shell intermediate states. In short, all hnite-volume effects arise from the sums over the 
two-particle loops shown explicitly in Eqs. (45) and ( [4^ . 


Einally we explain the functions A^^^'ik') (fc), focusing our wording on the former. A''^^^,{k') is the sum of 

all diagrams which contain exactly one insertion of a vertex factor, • • •), and two external legs. Additionally, 

the dehnition is restricted to only include diagrams which are two-particle irreducible in the (E,P) channel of the 
outgoing two-particle pair. The hnal restriction demands further explanation and is best summarized through the 
condition 


j(0) 


(0) 




(OlA(cr)A:-' 


(2/)A:-“"(0)|0) [SAk)]-^ = 


J L 


4(0) 

^L,a 


+ / ^A^ZAk')Sa'ik')ML,a'Ak',k). (49) 

1./ 


Note that the left-hand side simply defines the sum of every possible two-legged diagram with a single insertion of 


Af^}{k, • • •). Thus, this relation defines A^aAk') as the partial sum of diagrams in which the two-particle scatterings, 
shown explicitly in the second term, are not included. 

Here we have allowed for the possibility that A^^\,{k') and B^°)J(fc) may have non-negligible finite-volume depen¬ 
dence. This will occur whenever some of diagrams within these quantities contain on-shell intermediate states. For 
the purposes of this work, we restrict attention to operators A and for which no such diagrams occur. We thus set 
A^^\,{k') and B^j^'^^{k) equal to their infinite-volume counterparts and drop the L subscript from now on. We revisit 
the restrictions on interpolating operators in Section jlVj below, where we make specific choices for the operators in 
order to derive our main results. Since the diagrammatic representations of (fc') and will in general 

depend on the specific process in mind, in Figure we give examples of diagrams contributing to B^^'^^{k) for a 
particular reaction, N^* Ntt. 

We close this section my noting that spin has played a very minor role in the discussion. It is formally included 
everywhere above, and we stress in particular that that Eqs. (45l-(48) include implicit field indices. The functions 
A^^\, and 5^°)/ have one such index for each particle and the functions Ka'^a, Ba have two for each. The indices 


i(o) 


are thus contracted in the natural way, with viewed as a row vector, Ka'^a and Sa as matrices, and B^^^^ as a 

column. Similarly, the repeated channel indices are summed in all terms above. The new spin index structure and the 
modihed forms of the propagators are the main differences that arise here as compared to the case of scalar particles. 
We will show how to accommodate this change in the next subsection. 
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-C(P) F(P,L) TIHP) 


FIG. 2: As discussed in the text, the difference between the finite and infinite volu me two-particle loops is due to the cut of 
the loops, and can be written as the finite-volume matrix F{P, L), defined in Eq. (321, contracted with row and column vectors 
as shown. 


B. Finite-volume residue of two-particle loop 


Turning now to the second step listed above, in this subsection we identify the form of the finite-volume residue 
factor F{P,L). The factor generically appears in expressions of the form 


Na 




a=l 
No r 


L3 


a—1 k 


I ^Ca{P,k)Sa{k)TZliP,k) , 

l^Ca{P,k)Sa{k)nl{P,k) . 


" dk 


(50) 

(51) 


Here the functions Ca and TZj^ can stand for either the interpolating functions and or the kernels K. This 

notation is only used within this subsection, and we rely on context to distinguish these functions from the residue 
matrix and the Lagrangian. We now consider this finite-volume residue in three specific cases: two scalars, one scalar 
and one spin-half fermion, and two spin-half fermions. We then generalize the results for two arbitrary spin particles. 


1. Two scalar particles 


In the case of two scalar particles the finite-volume residue of the two-particle loop is given by 

Nr 




a—1 

Nr 


1 

Z3 


J ^ Ca{P,k)AUk)Al,{P - k)TZi{P,k), 




a=l 


L3 


>Ca(P,k:) 


2 uJal 2 uJa 2 {E - Ual - 0Ja2 + ie) 


-rKiPXa)- 


(52) 

(53) 


To reach the second equation here, we evaluated the integral over ki and discarded all contributions which are smooth 
functions of k. These give exponentially suppressed corrections, which we neglect throughout. We have also introduced 
new notation for the functions Ca and TZa, indicating that these only depend on total momentum P as well as the 
vector k which we have boosted to the CM frame. 

Next we use the identity proven in Ref. [22] which states that one can make the replacement 


where 


CaiPXaXiPXa) ^ CaTm.iP) 


Q.a 




,(P). 

a' ,1' ,m, ^ ' 


Ca{P,qlka) = Ca,l,mAP)^yimAK) , 
Tl\,{P.qlka) = 7^I,,;,,^;(P)V^F;^;(kI) ■ 


(54) 


(55) 

(56) 


We deduce 


Tj. = -C 


{P)FSm„a'l'ml (P, L)K',1 '(^) - 


(57) 
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where i® defined in Eq. (331 above. 

As mentioned above, for two-particles with momenta k and P — k subject to the constraint E = Wai -I- uJa 2 , the 
magnitude of a single particle’s momentum in the CM frame must be q*. Thus the key difference between the 
functions Ca{P, klk^) and Ca{P, g *k„ ) is that the latter has been projected on-shell. This projection emerges because 
the sum-integral difference in Eq. (|53[ ) is dominated by the two-particle pole. Equivalently, the difference between the 
left- and right-hand sides of Eq. (|54|) regulates the pole in Eq. (53). Thus the difference only generates exponentially 


suppressed corrections. Since we neglect these, the substitution is justified. 


2. One scalar and one spin-half fermion 


In the case of one scalar and one spin-half fermion the finite-volume residue takes the form 


ATc 




L3 


Pa{P,K) 


rriai 


k4 — iuJal 


2uJal‘2,UJa2{E — UJal — Wa2 + *e) 


PliPXa)- 


(58) 


This expression assumes that the propagating fermion has the quantum numbers of 'I'|0) (rather than drlO)). The 
distinction enters through the sign of ki = ioJai appearing in 
We now substitute the identity 


+ = E “aAi/2(k)Ua,Ai/,(k), 

Ai/ 2 =± 


(59) 


where MQ_Ai/ 2 (k) and Ma.Ai/a)^) are spinors used to interpolate states with definite spin. Here we are using helicity 
states, but this identity holds for spinors u and u, which interpolate spin states quantized along any fixed direction 
in the CM frame of the single spin-half particle. We deduce 


No 


-Al = - E 


1 

Z3 


^..A,,2(^.k:) 




/2 


2uJal2uJa2iE - UJal “ UJa2 + *e) 


Kx' 


(60) 


where 


CaMn (P^ k:) = Ca{P, k:) Ua,xa^) , (P, k:) ^ fia.A,/ 2 (k) 7^1(P, k:) . (61) 

Note that £a,Ai/ 2 (-P)ka) and (P, k*) are the definite spin projections of £a(P, k*) and Pjj(P, k*). For 

example consider the special case 


CaiP,K)= lim [ ’^'>y{0\TA{0)'i>{x)(j){y)\0)[iji:-\-mai][{P - kf +ml^] . (62) 

^4—,P4 —J 


Then it follows that 


CaMj.iP^A) = ( 0 |A( 0 )|P,P,k„s“,Ai/ 2 , in) . 


(63) 


The contracted spinor in Eq. (61) is precisely the required element to convert the correlator (evaluated in the LSZ 
limit) to a matrix element involving a physical two-particle state with definite helicity. We next decompose in spherical 
harmonics to reach 


P'a,Xi/2 (P’ 'laka) ~ P-a,l,Tni,Xi/2 (ko) ) 

^IcA(,,(^.9:k:)=7^Iq.,^;.A' ,(^)V4)ry;^,(k:). 


' 1/2 


(64) 

(65) 


As already discussed for general particles in Eqs. (18 1 and (HI) above, these may also be defined in a basis of definite 
total angular momentum 


Uj.mj(p)s Y1 {I Till , 1/2 ^1/2\J ,\i/2 i.^') ’ 

^^■1^1/2 

^(Im, V2A1/2IJ M)Ca,J,M,l[P) = P a,I,mi ,Ap/2 (P)- 


( 66 ) 


JM 


(67) 
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Applying the identity of Eq. (54) used for scalar particles above, we deduce 

(-P) L) ,l\m[ ,X\^ ' 

= -Pa,J,M,liP) FaJMl,a'J'M'l'{P,L) ^1'^M',1'’ 


( 68 ) 

(69) 


where 




{P,L) 


FaJMl,a'J'M'l'{P,L) = Ai/2 | JM) (i', I /2 Ai/2 | (P, L). 

iTii,Ai/2,mj 


(70) 

(71) 


These are just the definitions of Eqs. (31)-(331 above, applied to the special case of one scalar and one spin-half 
fermion. 

We stress here that the only new ingredient for spin-half particles relative to scalars is the on-shell Dirac matrix, 
which can be trivially accommodated using Eq. (59). The following points are crucial: 


It is justified to use the free form of the propagator in Eq. (58) only because of the sum-integral difference 
appearing in this equation. This allows one to replace fully dressed propagators with free propagators, because 
the difference between the two is a smooth function which thus gives only exponentially suppressed finite-volume 
corrections. 


• The observation that the on-shell Dirac matrix is equal to a sum of spinors, Eq. (591, is equivalent to noting 
that it is simply the identity matrix when written in spin space. This point is a clear physical necessity: The 
propagation of free-particle states must be independent of the azimuthal spin of those states. 


Using Eq. (59) also makes it clear that the specific value of the matrix in a particular basis of the gamma 
matrices is irrelevant. For example. Refs. [SOI [SI] first use symmetry arguments to set k • 7 to zero in Eq. ( |58| ) . 
This simplification requires that the remaining factors in the summand are invariant under k — > —k, and thus 
only applies in the case of P = 0. 


3. Two spin one-half fermions 

In the case of two spin-half fermions the finite-volume residue takes the form 


Na 




a—l 


1 ^ 


Pa,o.APXa) 

[-i^ + [- i{f -^) + ma2]^pl, 

p^ \ 


2uJal20Ja2{E - Wal - UJa2 + if) 


(72) 


Here we have explicitly shown the Dirac indices: a, a', (3, /?'. We do so to stress that the two matrices in the 
numerator act on different spaces. Substituting the decomposition in spinors, Eq. (59) above, we find 


Na 




1 

1 ? 


Ca,X.^,xAPXa) 




2 Wal2a;a2(£' - Wal - Wa2 + *e) 


<A' .A' (7^>k:) 

’ Si’ S 2 


where 


Pa,Xs^ .A.2 {P, k*) = £a,a,l3{P, k*) Ua,a,A,,(k) Ma./3.A,2 (? “ k) , 

.,A' ,A' (^>k*) = ^2a,a^A,^(k)Ma./3^A3,(P-k)7^^ ,^,(P,k*). 

’ Si ’ S 2 ^ ^ 


7^. 


In summary the helicity projection used above now appears twice for the case of two spin-half fermions. 
Decomposing in harmonics and applying the on-shell identity we find 

Pl £a,l,7ni,Xsi,XB2 Pal7niXsj^Xs2,a,'l'm'^XB^X'^J^P^ P) P'a' ,V ^m'^^X' ,A' ’ 

= —£a,J,ALl FaJMl,a'J’M'l'{P, L) ji M' l'{P) ) 


(73) 


(74) 

(75) 


(76) 

(77) 
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Cl(P) - C'oo(P) 



ui\ + ... = -A{P) — 


F-^{P,L)+MiP) 




FIG. 3: As discussed in Section [III C[ the difference between the finite and infinite volume correlation function is a geometric 
series in —M{P)F{P, L), where M{P) is the scattering amplitude and F{P, L) is the finite volume function defined in Eq. (321. 


where 


FajMis,a'J'M'i'S'{P, L) = 6ss' ^ {I "^ 1 , Sms\JM){I' m'l, S' ms\J'M'){P, L). (79) 


mi , 771 ^ ,ms 


Again these are the definitions of Eqs. (31)-(33) above, now applied to the special case of two spin-half fermions. 


4- General Spin 


Turning to general spin particles, we now make two key observations, true for all particle types, which allow us to 
reduce Eq. (50). Eirst we note that evaluating the integral and discarding terms which are exponentially suppressed 


generically gives an expression of the form 

Na 




1 

L3 


Ca{P,k) 


Ta{^) 


2uJal2uJa2iE - UJal “ + *e) 


PiiP,k), 


(80) 


where T is a matrix with index space of the two-propagators, i.e. two sets of field indices. Next we observe that 

/:,(P,fc)T,(k)pt(P,fc) nl^,^ yjp,k). (81) 


Here the factors of TZ and appearing on the right-hand side are redefinitions of the original C and TZ^ appearing 
on the left which couple to states with definite Agj. Further, these objects are normalized so that the infinite 
summations in the next section will give standard infinite-volume matrix elements and scattering amplitudes. As 
already mentioned above, we know that T must be proportional to the identity matrix for particles of any spin. This 
follows from the physical observation that for free-particle states the propagator cannot vary for different azimuthal 
components. 

Combining this result with the on-shell projection described by Eq. ( [^ , we deduce 

Fl = -C{P)F{P,L)TZ\P), (82) 


where C{P) is understood as a row-vector and TZ^{P) as a column vector with indices of either {Z} or {J}. F{P,L) 
is defined in Eqs. (31 )-(33l and must be used here with the same basis as £{P) and TZ'^{P). 


C. Summation of diagrams 


We now apply the key identity of the previous section 

N, 




a—l 


1 

L3 


J ^Ca{P,k)Sa{k)TZUP,k) = -C{P)F{P,L)TZ\P) 


(83) 


to the skeleton expansion given in Eqs. ( |45[ ) and (46). Each sum is replaced with an integral plus a residue which 
contains F. Reorganizing by number of F insertions, as shown in Figure]^ we deduce 


Cl{P) = Coo{P) + ^ A(P)[-E(P,L)][-M(P)F(P,L)]"Ht(p) ^ 


n —0 


(84) 
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and thus conclude Eq. ( [^ . 

As already mentioned in Footnote this result assumes that the correct scattering amplitude and matrix elements 
are given by summing the perturbative definitions to all orders. In particular the quantities Ap}(P), and 

A4{i}qp}(P) are reached by first defining the functions 

A{P,k) = A^°\P,k) + J-0^A^^\P,P)S{k')M{P,k\k), (85) 

Pt(p,fc') = i3(0)t(p,fc') + I^M{P,k',k)S{k)B(<^'>HP,k), (86) 

MiP,k\k) = -K{k\k)- J0^Kik\p)S{p)M{P,p,k), (87) 

which each have implicit indices a, s“, , S 2 j ^S 2 ■ Projecting the vectors k and k' on-shell reduces the coordinate 

^ ^ a. 

dependence to the unit vectors and . Finally, decomposing these in spherical harmonics gives A{q(P), P|j,j(P) 
and Alppp/}(P). 


IV. QUANTIZATION CONDITION AND MATRIX ELEMENT RELATIONS 


In this section we show how Eq. (351 can be used to derive a relation between finite-volume spectrum and scattering 


amplitudes, as well as relations between finite- and infinite-volume matrix elements. The main results of this section, 
and of the paper, are the matrix element relations, Eqs. (Ill I, (116), (119) and (123). 


Beginning with the relation between spectrum and scattering, we note that poles in C'i(P) are located at the 
energies of the finite-volume theory. More precisely, poles are located at P 4 = iPn,p,L where Pi,p,l, F2,p,l, • • • is the 
finite-volume spectrum at fixed P,L. Since no poles appear in C'oo(P), A{P), or B^{P) the spectrum is given by all 
energies for which the matrix 


P-i(P,T)+M(P) ’ ^ ^ 

has a divergent eigenvalue. Equivalently the spectrum is given by all zeroes of the function 

A(P,L) = det[p-i(P,P)-HM(P)]. (89) 

This is the fully general two-particle quantization condition which was first presented in Ref. [29) . 

We now turn to relating finite- and infinite-volume matrix elements. To accomplish this we consider the Fourier 
transformed correlator 


C'L(a:4 - 2/4,P) = dx dy e {0\TA{x)B\y)\0) 


= jdx j^dy ^ [(0|M(cr4, x)\E^, P, L)] ^ [(E„, P, L|S^(?/ 4 , y)|0) 


= j dx j dy 


-E^{xi-yi) 


( 0 |AI( 0 )|£;„, P, L)J ^ [(E„, P, l\bHo)\o) 


J L 




-Enixi-yi) 


(0|Ai(0)|i;„, p, L)\^ [(£;„, p, l\bHo)\o) 


(90) 

(91) 

(92) 

(93) 


In the second step we have assumed X 4 > 2/4 and have inserted a complete set of normalized finite-volume states. As 
a result of the Fourier transform, it is sufficient to insert only states with total momentum P. In the third step we 
have pulled out the time-evolution and translation operators. The time evolution operators give the exponential time 
dependence and the translation operators cancel the phase factors from the Fourier transforms. We thus conclude 
that the integrand does not depend on x or y, so that the integrals simply give factors of volume as shown in the 
final step. 
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FIG. 4: Shown are examples of diagrams that are present (a) as well as those that do not contribute (b) to for, for 


example, N^* —>■ Ntv. Equation (491 defines for generic systems. 


We now evaluate the same correlator in a different approach, by using Eq. (351 


Cl{x^ - 2/4, P) = L^ f , 


dPA 


,iP4{x4-y4) 


27r 


C^{P)-AiP) 


F-^{P,L)+M{P) 


bHp) 


(94) 

(95) 


= L3 

= E®”''"("'”"')^3(0|^(0)|£;„,P,{J},in) 7^p}.pq(^;„,P) (£;„,P,{J'},out|et(0)|0). (96) 


where 


P{j}.{j'}{En,P) = lim 

Pa^iEji 


— {iPa + En)-^ 


F-^{P,L)+M{P)\y^^y,^ ’ 


(97) 


is the residue of the matrix that appears between A{P) and B'''{P), evaluated at the nth two-particle energy. To go 
from Eq. (951 to Eq. (96) we assumed that xa > yA, allowing to close the contour in the upper half of the complex 


Pa plane. The only analytic structure encircled in the contour is the tower of finite-volume-spectrum poles along the 
positive imaginary axis. The integral thus reduces to a sum of residues at these poles. 

Equating Eqs. (93) and (96) we deduce 


(0|^(0)|£;„,P,L) (E„,P,L|6t(0)|0) 


J L L 


J L 


— (0|^(0)|i?„, P, {J}, in) (i?n, P)1 {En, p, {J'}, outlet (0)|0) . (98) 


From this result follows all matrix element relations presented in this work [Eqs. (Ill), (116), (119) and (123) 


below]. The equality relates matrix elements between finite-volume states and two-particle asymptotic states. We 
stress here that the result is only valid for E* = [E^ — below the lowest three or four-particle threshold. 

This restriction arrises because it is only for the energy poles below inelastic threshold that Eq. (35) is valid. Thus, 


although we know that the correlator in Eqs. (93) and (96) must equal an infinite series of decaying exponentials, we 
only quantitatively control the exponentials corresponding to two-particle states. We stress however that matching 
between (93) and (96) is unambiguous for each coefficient of the Euclidean-time-dependent exponentials. For this 


reason, as long as E* is below multi-particle threshold, then Eq. (98) includes all power-law finite-volume effects and 
only ignores exponentially suppressed corrections of the form e~™E 

Before deriving our main results from Eq. (98), we consider an alternative form of the relation. We begin by taking 
the ratio of the equation with a slightly modified version in which is replaced with . This gives 


(E„,P,L|^t(o)|o) 

L 

(0|^(0)|E;„,P,{J},in) 

E{J},{J'}{En,P) 

(,E„,P,{J'}, outlet (0)|0) 

■(£;„, P,L|St(0)|0)' 

L 

(0|^(0)|E„,P,{J},in) 

P{J},{J'}{En, P) 

(E„,P,{J'}, outlet (0)|0) 


Next we recall a result demonstrated in Ref. that the matrix TZ has only one-nonzero eigenvalue, and is thus 
equal to an outerproduct of vectors 


P{J},{J'}{En,P) 


cin cout 


( 100 ) 
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where 'S’| jj is understood as a column vector in space of open channels, and as a row vector. Substituting this 
into Eq. (991 gives 


{En,P,L\AH0)\0) 

^ (e;„,p,{j'}, outlet(0)|0) 

■(E„,P,L|Ht(0)|0)' 

(£;,,P,{J^},OUt|^t(0)|0) 

L 


;vt 

'^{J},{J'}(A'n,P) 

(E„,P,{J'}, outlet (0)|0) 

;vt 

^{J} 

'^{,/},{J'}(F^n,P) 

(E„,P,{J'}, outlet (0)|0) 


( 101 ) 

( 102 ) 


To reach the second equality we have multiplied the numerator and denominator by where is a 

completely arbitrary vector that can be freely chosen to make the relation as convenient as possible. The advantage 
of the second equality is that, although one can readily show that TZ can be written as an outerproduct, it can be 
complicated to determine the specific forms of f | "j and jq • 


A. One-to-two transition amplitudes with arbitrary spin 


In Ref. |27j . we demonstrated how to relate matrix elements of external currents to transition amplitudes by 
formally determining two-point and three-point functions and then taking appropriate ratios. In that earlier work we 
also discussed in great detail the dependence of the correlation functions on off-shell scattering amplitudes. In this 
section we circumvent a great deal of the discussion of our previous work and use Eqs. (98) and (1011 to generalize the 


matrix element relations to arbitrary spin. Beginning with (98), we note that one can formally choose the operators 
A and to satisfy 


B\x) = 
A{x) = 


■ -/ = :Ja(x) lim 

a/2Eo,pL3 ' Pi^^Eo 


P 4 —>ii5o 


+ M^] J {x + y), 

liin [P" -h M^] f d^ye-^^y^{x + y)jl{x). 


(103) 

(104) 


where Ja is an arbitrary local current and d> is the interpolating field for a single particle which is stable under the 
interactions governed by the Lagrangian, Eq. (36). 


H ere it is important to note that these specific choices for the operators satisfy the general for ms a ssumed to reach 
Eq. (35). To see this we must consider the endcap functions A^^\,(k') and i?^°|/(fc), defined in Eq. (49) above. For this 


discussion we actually define the endcap functions in terms of yj2EQ^plAA and ^J2EQ_pL'^B^ so that the prefactors 
in Eqs. (103) and (104) do not confuse the arguments. For the specific choice of Eq. (103), is equal to the 

sum of all one-to-two diagrams with the outgoing two-to-two finite volume scatterings amputated [see Figure for 
examples and a counterexample]. All contributing diagrams have no on-shell intermediate states and it thus follows 
that {k) is equal to its infinite-volume form, up to exponentially suppressed corrections. To see that this is the 
case, one must discount on-shell states both before and after the current insertion. Those appearing before are ruled 
out by the assumption that the incoming particle is stable. Those after the current must be two-particle states due 
to the restriction on total four-momentum. These are precisely the states that are amputated in the definition of 
Finally, for A^^\,(k') two-to-one diagrams must be considered. Applying the same arguments, we deduce 
that no on-shell states appear and thus that the finite-volume corrections to both objects are negligible. 

Returning to Eq. (103) and (104), the Fourier transform, propagator amputation and on-shell limit are the pre¬ 


scription required so that the operator picks out a one-particle state. In finite volume we find 

(E„, P, L|Hi(0)|0)] ^ = (K, P', L|:7 a(0)|Eo, P, L, 1), 
where the finite-volume one-particle state is normalized to unity 

(e;o,p,l,i|e;o,p,l,i) = 1. 

In infinite volume the result is 

1 


(E„,P,{J'},out|Ht(0)|0) = 


v/2Eo,pL3 


(E^, P, { J'}, out|,7A(0)|Eo.P, P, { J}i) 


(105) 

(106) 

(107) 































17 


where the infinite-volume states satisfy standard relativistic normalization 

{Eo^P-,P',{J}i\Eo,P,P,{J}i) = - P'), 


see 


also Eq. 0 abo ve. 


Substituting Eqs. ( |103[ ) and ( |104[ ) into Eq. ( |98[ ), we deduce 

|«,P',L|J^(0)|Eo,P,L,l)l = ^ ^ 


L3 




0,P 


n 


A.,{J} 




^out 


( 108 ) 


(109) 


where 


^X{j'}(^o.p, P; P') = (i?;, p', { J'}, out|jA(0)|Eo,p, P, {J}i). 
This can be trivially rewritten with all operators in momentum space 


IP', L\Ja{0, P - P')\Eo, P, L,l)\ = 




0,P 


n 


A,{J} 


7^{J},p},(£;;,p') 


njout 
A,{J'} ’ 


where 


[^XX'l (^o.p, P; K, P')] (2^)3<5XP - P' - Q) = (i?;, p', U'}, out\JAiO, Q)|Eo,p, P, {J}i), 
and the operator J'a is defined as the Eourier transform of the position space current, Ja, 

Ja{xo,Q) = J dxe~"^'^JA{x). 


( 110 ) 


( 111 ) 


( 112 ) 


(113) 


This is the most general possible Lellouch-Liischer relation for cubic volumes and two-particle states. We give an 
example of the utility of this result in Section [Vj and in Appendices and we discuss the free and narrow-width 
limits, respectively. 

Note next that Eq. (Ill) does not allow one to access the signs of transition amplitudes. Although the absolute 
sign of a given matrix element is not physically observable, the relative sign between two matrix elements is. With 
this in mind, we take Eq. (101) and substitute 


At(0) = ^ 

eXo) = 


1 


^^^^^^,(0) hin [P^ + M^] 
yy2EQ-pL^ P 4 —^^-Eo,p 

-=^^^,(0) hin [P^ + M^] 
yy2EQ-pL^ P 4 —>-^-Eo,p 



(114) 

rf4^g^Px^t(a.), 

(115) 


This gives 


(E;;,p^L|^AX0.P-P0l^o.P.^.l) 

(E;, P', LIJaXO, P - P')l^^o, P, A, 1) 


pout njout 
pout ^out 


;vt 

^{J} 


'TJOUt 

^Ai,{J'} 

;vt 

^{J} 


o/out 

^A2,{J'} 


(116) 

(117) 


where we used that in the ratio of finite-volume matrix elements, one can replace the position space current with a 
momentum space current, since these only differ by a factor of L^. 

Eq. (116) constrains the relative sign and may be more useful for other reasons in the analysis of a given physical 
system. In general, ^X*{j'} infinite vector over all open channels and all partial waves that mix in accordance 
with the symmetry of the system. If one ignores all but one entry, it is evident that there is a one-to-one mapping 
between the relative sign of finite-volume matrix elements and that between infinite-volume transition amplitudes. 
When including more than one partial wave or particle channel, Eqs. (Ill) and (116) allow one to simultaneously 
constrain the absolute value and relative sign of these transition amplitudes. This procedure is analogous to that 
which was implemented in a recent LQCD study of Ktt — Krj [5^ [53] . 
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B. Vacuum to two-particle transition amplitudes with arbitrary spin 

An even more straightforward result of Eq. ( [^ above is reached by setting 

= A(0) = 4(0), 

where, once again, Ja is a local current. One finds 


|(E„,P,L|J'^(0,-P 


= A V 


A,{J} 




V 


A,{J'} ’ 


where 


Similarly, substituting 


into Eq. (101) gives 


VX|jq(i?n,P) = (E„,P,{J'},out|^A(0)|0), 
[V°llj,y{E^,P)] {2TTfS^{P + Q) = {E^,P,{J'},out\JA{0,Cim. 

AH0) = JaA0), BH0) = JaM, 

{E^,P,L\JAA0,-Pm ^{jh 
(£;„,P,L|j7a,(0,-P)|0) £{% ’ 


7 ',, 

P{J}.{J'}iPri,P) 

■\;out 

7’n 

'^{J},{J'}(^^n,P) 

■\;out 

‘^^2.17'} 


(118) 

(119) 

( 120 ) 
( 121 ) 

( 122 ) 


(123) 


V. ONE NONTRIVIAL EXAMPLE: N + J^ {Nn, Nrj, Nrj’, T,K, AK) WITH |d| = 1 AND £ < 1 


Eq. (331 is a concise definition of the finite volume function When considering a specific example it is more 
convenient to introduce an alternative representation first used in Ref. |22| . 


na 


Klm^a'V^'SP^L) = 




I (k,)Ec.S (^a) 

l",m" 


, (124) 


where the functions are defined via 


cf^(fcf;L) = ^(^^y 'z,t[l;(fc*L/27r) 




= E 

rSfd 


\rnUr) 

(^2 _ ^2y 


(125) 


The sum is performed over Vd = {r G 
vector d = PL/27r, = 


1 + 


r = 7 ^(m— a£,d)}, m is an integer triplet, d is the normalized boost 
and 7 “^x = 7 “^X|| 


+ xj_, with 7 = E/E* and with xi i (xj_) denoting 


the X component that is parallel (perpendicular) to the total momentum, P. In Appendix]^ we demonstrate how to 
generalize this to describe systems in asymmetric volumes and systems with twisted boundary conditions. 

In Ref. [27j . we demonstrated how to utilize the main result presented here, Eq. ( |111[ ), for various systems with 
zero intrinsic spin. In order to illustrate the power of this generalization, we consider a simple scenario oi Ntt near 
threshold where contributions due to partial waves with £ > 2 can be ignored. Due the nonzero intrinsic spin, there 
are three partial waves we must consider 6ij = {<^Si/ 2 )'^A/ 2 : ^^ 3 / 2 }- When the system is at rest, parity is a good 
quantum number and the S-wave and P-waves do not mix. When the system has nonzero total momentum, parity is 
no longer a good quantum number and consequently S-wave and P-waves can in fact mix. This mixing is explicitly 


illustrated in the finite-volume functions appearing in the quantization condition, Eq. (Ill), and the residue matrix, 
Eq. (1^. 
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For example, when the system has a boost vector |d| = 1, its symmetry group is Dic 4 . Following the notation used 
in Ref. [54], the irreps that couple to half-integer spin systems are the Ei = S'i/ 2 ©Pi/ 2 ©F 3 / 2 ® ■' • and E 3 = P 3 / 2 © - ■ •. 
The finite-volume functions and scattering matrices corresponding to these irreps are 


Dic4 El 


Fw, = 


SnE* 


A4ei = 


Dic 4 E 3 


8nE* 
q* 

* 

F^ = ^ 
8nE* 


cot ((iqq + z 

cot 

©3 

01 cot (j)% 

\ 

cot 05*0 

V3 

cot + i 

-0§COt (j)^Q 


cot 

-01 cot (j)% 

cot </>^0 + ^ 

+v 

[cot Ss, - i] 

-1 0 

0 

\ 

0 

[coti^Pj - 

7]-* 0 


0 

b 

[cot(5p3 — 
2 

rV 


d cot (/)fo 
cot (j)oo 


where the pseudophases, ^^e defined via 


75 


47r 


+ i 


Me. = 


8ttE* 


1 


q* cot dpj — V 


9a, n cot L). 

9a.« 


(126) 

(127) 

(128) 

(129) 


Determining these matrices as functions of E and L and substituting into Eq. (971 gives a three-by-three matrix, 
"TT-Ei (E'Ei,n, P), and a single value, 77,E3(pE3,n, P), for each energy level in the finite-volume spectra of the indicated 
irreps. The single value, TZ-e^^Ee^^u, P), gives a Lellouch-Liischer like proportionality, relating the finite-volume matrix 
element of the nth state to the transition amplitude 


I (p;, p', L, E 31 77^(0, p - p')|i7o.p, p, m)\ = 


\/‘^Fqi 


7^E3(E;,P')| \HZp,\, 


where 


[HTk (i^o.p, P; E'„, P')] (2^)353(P - P' - Q) = (p;, P', N, tt, Pa, out|J)i(0, Q)|Po,p, P, N) . 
2 

In this case one can easily simplify the form of TZe^ 


7?.e3(E„,P) = 


d 


1 -1 


{F-\P,L)+MiP)) 


dE 

q 


E^Er, 


d 


1 -1 


M^iP)^{FiP,L)+M-\P)) 


E=E„ 


8ttE* 


sin^( 5 p 3 e * I cot (/(qq — 


cot (j)' 


20 


q „“2ii5p3,2 


8'nE' 


JL 

We 


dE 

(j^Es + ^P| ) 


75 


cot Sp, 


E=E„ 


J E=E„ 


where we have introduced 


(H (H cot 09(1 

cot (j)^^ = cot doo - 


75 


(130) 

(131) 

(132) 

(133) 

(134) 

(135) 


Note that the phase appearing in 77 .E 3 is exactly that needed to cancel the phases from 77“ and 77°“* so that the 
right-hand side is pure real. 

For the three-by-three matrix, TZei, there is no straightforward reduction of the general relation, Eq. (111). In this 
case the result relates the finite volume matrix elements of states that transform in Ei to linear combinations of three 
transition amplitudes 


[77X7(^o.p,P;K,p') 


(27r)3d3(P-P'-Q) = 


/(p;, p', tv, , out|^A(0, Q)|Po,p, p, iv)\ 
(p;, P', TV, TT, Pi , out| J^(0, Q) iPo.P, p, TV) 
V(p;, P', TV, TT, P|, out|J^(0, Q)|Po,P, p, TV) J 


(136) 
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For physical or nearly physical quark masses, the three-particle {Nttit) threshold resides close to the Ntt threshold. 
For unphysically heavy quark masses, the three-particle threshold is pushed up and other two-particle thresholds 
{Nt], Nr]', YjK, AK) approach the Ntt threshold. In the SU(3) flavor limit all of these thresholds overlap. As we have 
discussed, incorporating more than one two-body channel amounts to upgrading the angular-momentum matrices to 
matrices in the product space of flavor and angular momentum. 


VI. CONCLUSION 


In this work we present a relation between finite-volume matrix elements and 1 —>■ 2 as well as 0 —> 2 transition 
amplitudes in the presence of an external current. The result is exact up to exponentially suppressed volume corrections 
and is nonpertubative in the strong dynamics. The result presented here is the most general of its kind, and holds 
when the individual particles have arbitrary spin and the final state is composed of any number of strongly coupled 
two-body channels. Furthermore, the result is independent of the nature of the external current. It may therefore be 
implemented for determining electroweak transition amplitudes as well as BSM processes. 

To study transition amplitudes involving two particles in the initial or final state, one must first constrain the 
scattering amplitude. This can be done by determining the finite-volume spectrum and then applying Eq. (0/@ 
as has been done, for example, in Refs. [SailalES]. Once the phase shift and mixing angles have been parametrized, 
one may proceed to take derivatives of these to determine the residue function TZ, Eq. , which allows one to 

relate finite-volume matrix elements with transition amplitudes. In general, in order to determine transition processes 
involving resonances, it is necessary to determine not just the matrix element of the ground state but also excited 
states. Recently, Ref. [2] has demonstrated an efficient way to determine excited state matrix-elements using the 
distillation framework. Having constrained both the residue function and the finite volume matrix elements for 1 —> 2 
and 0 —» 2 processes, one may proceed to access the corresponding transition amplitudes using Eqs. d^/dml) and 
Q/(119), respectively. 
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Appendix A: Free limit 


To give further insight to Eq. (Ill), here we first consider the result near two-particle production threshold and 
then in the case where the two outgoing particles are non-interacting. In the former case, we set all phase shifts to 
zero with the exception of the S-wave. In this limit, the ratio between the transition amplitude and finite volume 
matrix element simplifies to 




16ttEo^pE': + 0 


|(£;;,P',L|J(0,P-P')|Ao,p,P,L,l)|2 CMAA/ q- 


BE' 


(Al) 


E'=E' 


where, for clarity, we have explicitly included the normalization factors of the initial (A/i) and final (A//) finite-volume 
states. Next, using Eq. (129), one obtains a simple result in the limit in which the final two particles do not interact, 




2i?o,P 1 ^ 


2AoP E'^ 


|(f;;,p',l|j(o,p-p')|Ao,p,p,l,i)| 2 ^ MAf/ e, 




(A2) 


where is the degeneracy of the nth state, equivalently the number of vectors n S such that = n. In general, 
one expects TZ/TZfree ~ 0{\) for weakly interacting systems. In arriving at this result, we have used the fact that the 
phase shift, and thus also its first derivative, vanishes, as well as the identity 84)^^ = cos^ 9tan(/)QQ. 
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FIG. 5: (a) In the vicinity of a resonance, the two-body scattering amplitude can be obtained by mediating the two-body 
systems using a fully dressed auxiliary field, (b) Using this same auxiliary field, one can also obtain an expression for the 1 —> 2 
transition amplitude in the presence of an external field, Eq. (B4|. 


Appendix B: Narrow-width approximation 


Another interesting limit, which was previously considered in Ref. |2f)j . is the narrow width approximation. Here the 
exact nature of the resonance is not relevant, and thus we leave the partial wave, £, unspecified. In the narrow-width 
limit, the two-body resonance approaches a bound state. A resonance corresponds to a pole in the complex plane, 
with the imaginary part of the pole location proportional to the resonance width. If we send the width towards zero, 
then the pole approaches the real axis and the state becomes bound. In order to consider this limit, we investigate 
the behavior of the scattering amplitude and transition amplitudes near a resonance. We begin by considering the 
two-particle scattering amplitude at energies near the resonance pole. Using a Breit- Wigner inspired parametrization, 
we write the scattering amplitude as [see Fig. |5(a)| 


Mi = 


GUniE*) 


- E*^ - iE*T{E*) 


(Bl) 


where 


is a generic function that parametrizes the coupling of the two-particle state to the intermediate 
resonance, vriR is the real part of the resonance pole location and r(£l*) is its energy dependent decay width. From 
Eq. (29) above, we know that the scattering amplitude must satisfy 


Mt = 


SttE* 


1 


^q* cot Si — i 

Equating these two expressions we reach the following relations for G 2 ^r, F and the scattering phase shift, 

G 2 ^r{E ) - , tan5(E ) - • 


(B2) 


(B3) 


Similarly, we define the transition amplitude in the vicinity of the resonance as [see Fig. |5(b)] 

G2^RiE*) 


ni = 


m\ - E*^ - iE*T{E*) 


= E^^^{E*,Q^) 


l87rE*'^r{E*) 


- E*"^ - iE*T{E*) 


iq* 


(B4) 


where ,Q'^) is a smooth function of both E* and = {P — P'Y. Note that although the numerator of the 

scattering amplitude in the presence of a resonance is proportional to r(£’*), the transition amplitude is proportional 
to 

Now we assume we are in the vicinity of a single resonance state. Using the parametrization of the phase shift given 
above and defining Fj^ = r(mij), one finds the derivative of the phase shift as it approaches 7r/2 to be 




E*—mR 


= — sin^ (5^^ cot 5(E*) 
oE 


= 2 

E-=mn ^RITT-r 


[1 - 1 - 0{TR/mR)] , 


(B5) 


where Er = y + rn^R- As the width goes to zero the derivative of the phase shift divergences as expected. Ignoring 
the contribution from the derivative of the pseudophases, which are finite away from free-particle poles, one obtains 
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the following result for the right hand side of Eq. (A1) near a narrow-width resonance 

2i?o,p IQiiEfi 


_ _ _ 

|(£;,i,P',L|J(0,P-P')|Eo,P,P,L,l)|2 


[1 + 0{VRlmR)] , 


(B6) 


From Eq. (B4) we find that near the resonance mass, the transition amplitude divergences inversely proportional 
to the square root of the resonance width 


ni = Fl^^{mR,Q") 


Stt 






(B7) 


In this limit, one finds that the finite-volume matrix element is equal to the infinite volume 1 
amplitude up to the standard normalization of the states. 


R transition 


\{Er, P', L|J( 0 , P - P')| Ao.p, P, L, 1 ) I" = \F}^^{mR, Q^)\ 


/02m2 ■N'iNf 


‘2Eqp2Er 




(B8) 


It is important to emphasize that this approximation only holds when the resonance is very narrow (Ep/mi^ ^1) and 
the energy level determined corresponds to the resonances mass, up to small corrections that scale with the width. To 
reliably asses the validity of this approximation, one must first determine the phase shift as a function of the energy 
using Eq. (891 as done, for example, in Ref. [55] for the p-resonance. Since a zero width resonance is equivalent to a 


two-body bound state, this discussion applies also for bound states, when the width is exactly zero. 


Appendix C: Asymmetric boxes with twisted boundary conditions 


Although the vast majority of present day calculations are performed in symmetric volumes with periodic boundary 
conditions |56j , the formalism presented here can be generalized to volumes that are arbitrary rectangular prisms with 
twisted boundary conditions. Following the notation used in Refs. nzmi], we define <pa,i and 4 >a ,2 to be the three- 
dimensional phases of the first and second particle, respectively, in the ath channel. We also define L to be the 
spatial extent of the z-axis and rji to be defined such that Lj, = r]xL and Ly = rjyL. Furthermore, let the vector 


X = {Xx/VxtXy/VyjXz)- With this, we can express the generalization of the functions in Eq. (1251 


2 „1 _ Ylmjr) 

^Im 1*1-^ Elx,ny\ — 2-^ (■j.2 _ 3,2)8 ’ 


(Cl) 

(C2) 




4>i,4>2'> 


where V 




= |r G I r = 7 ^(m- Oad + with m a triplet integer, = -(oa - + 0a, 2 ) + 


|(0a,i — 0a, 2 ), d = PA/27r, and all other quantities defined after Eq. (125|. 


By replacing the factors of F in all expressions with this generalization, and additionally replacing all factors of 
with rj^rjyL^, one reaches extensions of all results that are valid for arbitrary rectangular prisms with arbitrary 
twist. As discussed in Ref. in the presence of general twisted boundary conditions, the symmetry of the system 
is further reduced. Depending on the twists chosen, one may need to rederive the E-functions appearing in Eqs. (126) 
and (128). We point the reader to Refs. [271 mi [57] for further details. 
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